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Abstract 

In this survey we explore implementations of re¬ 
cently proposed distortion synthesis techniques, 
namely the Feedback Amplitude Modulation (1st 
and 2nd order cases) and Allpass filter coefficient 
modulation. These techniques are based on Periodic 
Linear Time-Varying systems, in which we operate 
by finding a suitable modulation function to obtain 
a desired spectrum. In order to illustrate this survey 
and encourage exploration of these new techniques 
we present examples in the Csound language. 
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1 Introduction 

Knowledge of a good set of synthesis techniques 
is a must for electronic / electroacoustic / com¬ 
puter musicians, sound designers and synthe¬ 
siser builders. Audio effects plugins developers 
also benefit from dealing with synthesis tech¬ 
niques that can also be used in an audio effect 
context, such as the ones that will be presented. 

In a previous survey, classic and recent dis¬ 
tortion synthesis techniques along with their 
implementations were presented at the LAC 
[Lazzarini, 2009]. Waveshaping [LeBrun, 1979], 
phaseshaping [Ishibashi, 1987], and summation 
formulae [Moorer, 1976], some of the classic ap¬ 
proaches to distortion synthesis, were analysed 
mathematically and demonstrated in Csound 
code. More recent or unusual approaches like 
asymmetrical FM [Palamin et al., 1988], Phase 
Aligned Formants [Puckette, 1995], and Modi¬ 
fied FM synthesis [Lazzarini and Timoney, 2010] 
were also exposed. Other, more common dis¬ 
tortion techniques, are also widely presented in 
computer music textbooks [Dodge and Jerse, 
1997] [Moore, 1990] [Puckette, 2007]. 

In this paper, we address some recent tech¬ 
niques which extend the distortion synthesis 
family, namely techniques based on Periodic 


Linear Time-Varying (PLTV) systems. With 
these new approaches we operate by modulat¬ 
ing the coefficients of a filter, and it was shown 
[Pekonen, 2008] that it results in a kind of 
dynamic version of phase distortion [Ishibashi, 
1987]. A thorough study of time-varying sys¬ 
tems applied to these techniques [Cherniakov, 
2003] [Timoney et ah, 2014] discloses the ap¬ 
propriate tools for properly understanding and 
using these new synthesis processes. These sys¬ 
tems are intended to operate differently than 
audio effects such as time-varying delay lines 
used in flanging [Zolzer, 2011] and allpass filters 
used to obtain variable fractional delays [Peko¬ 
nen et al., 2010]. Additionally, they do not re¬ 
place another type of variation within musical 
systems where users manually change (e.g. with 
a knob) parameter values to achieve a sonic ef¬ 
fect, such as mentioned in [Wishnick, 2014], 

Our motivation for presenting this survey 
comes from the fact that most audio program¬ 
mers’ work is supported by classic Linear Time- 
Invariant (LTI) systems theory [Oppenheim and 
Schafer, 1975], but Linear Time Varying (LTV) 
systems theory is less well covered in the liter¬ 
ature [Huang and Aggarwal, 1982]. Also, time- 
varying tools were indeed considered a long time 
ago [Layzer, 1971] [Risset, 1969] for audio appli¬ 
cations within our context, but only recently are 
being tackled in a more comprehensive fashion. 
A good understanding of LTV theory can bring 
us new kinds of synthesis/effects techniques and 
different ways for implementing established ones 
(which can bring nice variations or reduce com¬ 
putational costs). 

We will dedicate each of the next sections to 
one technique and then conclude the text. The 
aim of this paper is not to delve into the the¬ 
ory of LTV systems; our intention, instead, is 
to present another look at some techniques and 
also provide code 1 for their implementation. 


1 http://www.ime.usp.br/~ag/dl/lacl5-code.zip 



2 Allpass coefficient modulation 

An allpass filter contains poles and zeros at re¬ 
ciprocal distances from the origin, so their ef¬ 
fects on all frequencies are balanced [Moore, 
1990]. A typical 1st order invariant allpass 
would have a fixed coefficient, but a time- 
varying one has the transfer function 


In order to get an approximation of a sawtooth 
from a cosine, we must read the rising part of 
the cosine in a shorter time, and take more time 
to read the decaying portion. The function we 
must add to the phase generator, in this case, 
is drawn on the upper panel of Figure 1 and is 
given by [Lazzarini et ah, 2009b] 




— m(t) + z 1 
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with m(t) as the modulating function which 
drives the coefficient. Using this filter, a time- 
varying phase distortion given by [Timoney et 
al., 2009] [Laakso et ah, 1996] 
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is introduced in the signal. If we know how 
we want the phase to be distorted - in other 
words, if we know (j>(uj, t ) - we can use the ap¬ 
proximation tan (x) ~ x [Timoney et ah, 2009] 
and determine the modulating signal as 


m(t) 


_ -Q(AU) + uj) _ 

2 sin (cu) — t) + uj) cos (uj) 


( 3 ) 


The allpass filtering described can be imple¬ 
mented with the difference equation [Lazzarini 
et al., 2009b] 


y(n) = x(n — 1) — m(n)(x(n) — y(n — 1)). (4) 


The condition |m(n)| < l,Vn, assures stability 
[Cherniakov, 2003] and a DC offset 


DC(n) 


1 — m(n ) 
1 + m(n ) 


( 5 ) 


is introduced in the signal [Pekonen, 2008]. 


2.1 Classic phase distortion emulation 

In [Pekonen, 2008] [Lazzarini et al., 2009b] [Ti¬ 
money et al., 2014], we find as an example for 
this technique the emulation of the classic phase 
distortion [Ishibashi, 1987]. The phase distor¬ 
tion algorithm consists of reading a cosine table 
in an unusual way. Instead of getting indexes 
from a regular phase generator (which goes from 
0 to 1 in a period related to the chosen funda¬ 
mental frequency) multiplied by the table size, 
we add another function to the phase genera¬ 
tor values and then read the cosine wavetable. 


(I — d )f, x< 

x>d, 


( 6 ) 


where d is how long it takes to read from the 
start of the cosine table up to its maximum, 
and the smaller it is the more abrupt the rising 
ramp will be and the more distortion we will 
obtain. 

If we want to implement g(x) using a mod¬ 
ulated allpass filter we will get better results if 
i j>(uj,t) goes from — uj to — n [Lazzarini et al., 
2009b], so we make 
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The resulting modulation function is shown on 
the lower panel of Figure 1. On Figures 2 
and 3 we can see the waveforms and spectra 
of the sawtooths generated with both the origi¬ 
nal technique and the modulated allpass based 
one. Notice that the missing components of the 
classic phase distortion technique are actually 
present in the spectrum of the modulated all¬ 
pass output. We can also see and hear that 
the latter spectrum is richer. The code for im¬ 
plementing this example is given in Listing 1, 
where we can see two instruments; the first one 
works as a synthesis instrument, distorting the 
sinusoid and producing richer spectra; the sec¬ 
ond one applies the same technique as an au¬ 
dio effect, distorting the sound of a pre-recorded 
flute (any sound file can be used, or even a mi¬ 
crophone input). 


Listing 1: Classic phase distortion emulation as 
synthesis (instrument 1) and technique used as 
an audio effect (instrument 2) 

1 <CsoundSynthesizer> 

2 

3 <CsOptions> 

4 —o dac 

5 </CsOptions> 

6 

7 <CsInstruments> 

8 

9 0dbfs=l 
10 





Figure 1: Phase distortion function (up¬ 
per panel) and coefficient modulation function 
(lower panel). Source: [Timoney et al., 2014] 



Figure 2: Waveforms generated with the classic 
technique (solid fine) and the modulated allpass 
(dashed line). Source: [Timoney et al., 2014] 


JIllll 


1 1 1 1 1 11 1 111 . 


Figure 3: Spectra obtained with the classic tech¬ 
nique (upper panel) and with modulated allpass 
(lower panel). Source: [Timoney et al., 2014] 


11 /* audio —rate coeff allpass */ 

12 opcode Allpass , a , aa 


13 

adel init 0 


14 

setksmps 1 


15 

audio,acoef xin 


16 

aw = audio + acoe 

d* adel 

17 

ay = —aw*acoef + 

adel 

18 

adel = aw 


19 

xout ay 


20 

endop 


21 



22 



23 

/* PD function */ 


24 

/* inflection point 

*/ 

25 

gip = 0.1 


26 



27 

/* ftgen producing 

line S' 

28 

ipdfun ftgen 1,0,16384,7,1 


,1,16384*(1 — gip 

) ,o 

29 



30 

/* 


31 

instr 1: 


32 

PD synthesis 


33 

using sine wave in 

put 

34 

*/ 


35 



36 

instr 1 


37 

ifr = p5 


38 

iamp = p4 


39 



40 

/* regular phase 

genera 


*/ 


41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 
61 
62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 


aph phasor ifr 

/* phase distortion signal */ 
apd tab lei aph, 1, 1, 0, 1 

/* scaling of pd signal (eq.7)*/ 
iw = 2*$M_PI* ifr /sr /*omega value*/ 
ia = (1 — 2*gip)*$M_PI 
apd = apd*ia 

apd = apd*(ia — iw)/ia — ia — iw 

/* coefficient modulation function 
obtained from phase distortion 
signal (equation 3) */ 

/* envelope for modulation */ 
kmod linseg 0,1 ,1 ,p3 — 2,1,1,0 
amod =— kmod* (apd + iw)/(2*sin(iw) 

— (apd+iw)* cos(iw)) 

/* sine input */ 

asin tablei aph,— 1 ,1 ,0,1 

/* allpass */ 

asig Allpass asin ,amod 

/* envelope */ 

aout linenr asig * iamp, 0.01, 0.1, 
0.01 


outs aout 
endin 


aout 


/* 

inst r 2: 

PD adaptive synthesis 

using a monophonic instr input 











































74 */ 

75 

76 instr 2 

77 iamp = p4 

78 /* input signal */ 

79 afl diskin2 ” flutec3 . wav” , 1 ,0,1 

80 

81 /* pitch tracking */ 

82 kfr ,kamp ptrack afl ,2048 

83 kfr port kfr ,0.01 

84 

85 /* master phase signal */ 

86 aph phasor kfr 

87 /* phase distortion signal */ 

88 apd tab lei aph,1,1,0,1 

89 

90 /* scaling of pd signal */ 

91 kw = 2*$M_PI* kfr / sr 

92 ia = (1 — 2*gip)*$M_PI 

93 

94 apd = apd*ia 

95 apd = apd*(ia — kw)/ia — ia — kw 

96 

97 /* envelope for modulation */ 

98 kmod linseg 0,1,1,p3 — 2,1,1,0 

99 

100 amod = —kmod*(apd + kw) /(2* sin (kw) 

— (apd+kw)* cos(kw)) 

101 

102 asig Allpass afl , amod 

103 aout linenr asig * iamp, 0.01, 0.1, 

0.01 

104 

105 outs aout, aout 

106 endin 

107 

108 </CsInstruments> 

109 

110 

111 <CsScore> 

112 

113 /* uncomment lines to 

114 run instruments */ 

115 i 1 0 10 0.5 440 

116 ; i 2 0 10 0.5 

117 

118 

119 </CsScore> 

120 
121 

122 </CsoundSynthesizer> 


2.2 Finding a distortion function 

Now we present a new example with the allpass 
coefficient modulation technique to bring more 
insight. In this example, instead of choosing a 
desired waveform and then finding out how to 
modulate the allpass coefficient, we choose an 
arbitrary signal to distort the phase. Keeping 
in mind that we should generate a signal within 
the appropriate range, that is [-1,1], any signal 
a priori can be considered for the process. 

So we embarked on an experiment that in¬ 


volved the summing of partials in order to find 
a function to distort the phase of an input sinu¬ 
soid and thus get a more musically interesting 
timbre. The expression we used to create our 
example is given by 


yin) = 0.4 cos (/ 0 ) + 0.4 cos (2f 0 - 
+ 0.35cos (3/o + ^) + °- 3cos ^ 4 /o + -y 

( 8 ) 

In order to shift the output of Equation 8 to the 
appropriate range we use the scaling 

*<») = ( 9 ) 

Listing 2 presents Csound code for implement¬ 
ing this new example. The upper panel of 
Figure 4 shows the distortion function for this 
example, while the lower panel shows the all¬ 
pass modulation function resulting after apply¬ 
ing Equation 3. Figure 5 shows the waveform 
and spectrum obtained. Albeit this was a naive 
example, it demonstrates in an intuitive way 
how a phase function could be easily generated 
using additive synthesis and then used to drive 
the allpass filter coefficient after applying Equa¬ 
tions 9 and then 3, highlighting another musical 
possibility for using the allpass filter. 


Phase Distortion Function 



Samples Index 


Allpass Coefficient Modulation Function 



Figure 4: Phase distortion (upper panel) and 
resultant modulation (lower panel) functions. 


Listing 2: Implementation of an arbitrary func¬ 
tion as a phase distorter. 

1 <CsoundSynthesizer> 

2 

3 <CsOptions> 





Modulated Allpass Filter Output 

1 . 51 - 1 - 1 - 1 - 


45 

46 outs aout, aout 

47 endin 

48 

49 </CsInstruments> 

50 

51 <CsScore> 

52 i 1 0 10 0.25 440 

53 </CsScore> 

54 

55 </CsoundSynthesizer> 




»ii 11111111 ii i ii 111 ii 111 ii 11 ii ii ii ii 11 ii i ii ii 11 ii i 
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Frequency (Hertz) x ^ g 4 

Figure 5: Waveform and spectrum obtained 
with arbitrary function 

4 —o dac 

5 </CsOptions> 

6 

7 <CsInstruments> 

8 0dbfs=l 

9 

10 /* audio —rate coeff allpass */ 

11 opcode Allpass, a, aa 

12 adel init 0 

13 setksmps 1 

14 audio,acoef xin 

15 aw = audio + acoef*adel 

16 ay = —aw* acoef + adel 

17 adel = aw 

18 xout ay 

19 endop 

20 

21 /* PD function */ 

22 ipdfun ftgen 1, 0, 16384, 9, 1, 0.4, 

-90, 2, 0.4, -150, 3, 0.35, 
180/7-90, 4, 0.3, 4*180/3-90 

23 

24 instr 1 

25 

26 ifr = p5 

27 iamp = p4 

28 

29 iw = 2*$M_PI* i fr/sr ; omega 

30 

31 aph phasor ifr ;regular phase 

32 /* phase distortion signal */ 

33 apd tab lei aph, 1, 1, 0, 1 

34 

35 /* coef mod signal (eq.3)*/ 

36 apd = —0.5*$M_PI*(apd + 1)/2 ; scaling 

37 amod = — l*(apd + iw)/(2*sin(iw) — ( 

apd+iw)*cos (iw) ) 

38 

39 /* sine input */ 

40 asin tablei aph,— 1,1,0,1 

41 /* allpass */ 

42 asig Allpass asin,amod 

43 /* envelope */ 

44 aout linenr asig *iamp , 0.0 1 , 0.1 ,0.0 1 


3 Feedback amplitude modulation 

The Feedback Amplitude Modulation was men¬ 
tioned by Layzer (1971) and described and im¬ 
plemented by Risset (1969) in his catalogue ex¬ 
ample #510, but a rigorous mathematical anal¬ 
ysis of the technique was lacking. Since 2009 
there was a renewed research interest in ex¬ 
ploiting its musical possibilities [Lazzarini et al., 
2009a], [Kleimola et al., 2011], [Lazzarini et al., 
2011 ], 

First of all we will review its 1st order case. 
The basic idea is to modulate the amplitude of 
an oscillator using its previous output, as in 

y (n) = cos(w 0 n)[l + y(n- 1)], (10) 

with ujq = 2irfo and the initial condition y(n) = 
0 for n < 0. 

A first analysis [Kleimola et al., 2011] is made 
expanding Equation 10 as 

y(n) = cos (u 0 n)+ 

cos (wo n) cos ( u)o[n — 1])+ 
cos (won) cos (wo[n — 1]) cos (wo[rr — 2])+ 

... ( 11 ) 


oo k 

y{n) = X! IT c °s(w°[n-m]), (12) 

k =0 m =0 

showing that the resultant signal is composed by 
harmonics of the fundamental fo . A more inter¬ 
esting case is when we can control the amount 
of feedback in the system, so we introduce the 
feedback parameter (5 and Equation 10 becomes 

y(n) = cos (£n 0 n)[l + /3y(n — 1)]. (13) 

The feedback parameter can be interpreted to 
be similar to modulation index of conventional 
FM synthesis, and its influence on FBAM’s 





spectral evolution is shown in Figure 6, bor¬ 
rowed from [Kleimola et al., 2011]. 
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Figure 6: (3 influence on FBAM spectrum. 

Source: [Kleimola et al., 2011] 

The interpretation of the system as a time- 
varying filter enables a better analysis. The 
structure for the 1st order case is 

y(n) = x(n) + a(n)y(n - 1), (14) 

with 


x(n) = a(n) = cos (uon), (15) 

so we have a coefficient modulated HR filter. 
As in the previous section, the modulation will 
create the dynamic phase distortion, generating 
new partials, and in this case the filter is not 
allpass, so it has a time-varying non-flat magni¬ 
tude response. 

Applying the equations for stability analysis 
by Cherniakov (2003) a condition for stability 
of this system was derived as 


N 

pU cos (cuo rn) 

m= 1 


< i. 


(16) 


An expression for the maximum value of beta 
was proposed by Kleimola et al. (2011) as 


/3 m ax ~ 1-9986 - 0.00003532(/ 0 - 27.5), (17) 

and they also showed that as the value of f3 is in¬ 
creased, the presence of significant components 
in the output due to aliasing can be observed 
before the system becomes unstable. Listing 3 
presents code for the 1st order FBAM imple¬ 
mentation. 


Listing 3: 1st order FBAM 


1 

<CsoundSynthesizer> 


2 



3 

<CsOptions> 


4 

—o dac 


5 

</CsOptions> 


6 



7 

<CsInstruments> 


8 



9 

0 dbfs= 1 


10 



11 

opcode FBAM, a, kkki 


12 

; set vector size to 1 sample 

13 

setksmps 1 


14 

ay init 0 ; y[0] =0 


15 

ka, kf, kb , ifn xin 


16 

/* current sample + 


17 

weighted previous sample 

*/ 

18 

ay oscili ka + kb*ay, kf , 

ifn 

19 

xout ay 


20 

endop 


21 



22 

instr 1 


23 

; amp, freq , beta 


24 

al FBAM 0.5, 440, 0.7, 1 


25 

out al 


26 

endin 


27 



28 

</CsInstruments> 


29 



30 

<CsScore> 


31 

fl 0 16384 10 1 


32 

il 0 5 


33 

</CsScore> 


34 



35 

</CsoundSynthesizer> 



The 2nd order FBAM is obtained using two 
previous outputs in the modulation, each with 
its own feedback parameter. The system equa¬ 
tion is then given by 


y{ri) = cos (w 0 n)[ 1 + /3iy(n - 1) + f3 2 y(n - 2)]. 

(18) 

It was shown [Lazzarini et al., 2011] that with 
this system we can get a narrower pulse, and 
thus a richer spectral output for the FBAM sys¬ 
tem, as shown in Figure 7, borrowed from [Laz¬ 
zarini et al., 2011]. Code for implementing the 
2nd order FBAM is presented in Listing 4. 

Listing 4: 2nd order FBAM 

1 <CsoundSynthesizer> 

2 

3 <CsOptions> 

4 —o dac 

5 </CsOptions> 

6 

7 <CsInstruments> 

8 ksmps = 10 

9 Odbfs = 1 
10 nchnls = 2 


































Frequency (kHz) 


Figure 7: Comparison between 1st (in dots) and 
2nd order (solid line) FBAM for 500 Hz funda¬ 
mental. Source: [Lazzarini et al., 2011] 

11 

12 opcode Fbam2 ,a,akk 

13 setksmps 1 

14 asigml init 0 

15 asigm2 init 0 

16 ain,kbl,kb2 xin 

17 asig = ain*(l + kbl*asigml + kb2* 

asigm2) 

18 asigm2 = asigml 

19 asigml = asig 

20 xout asig 

21 endop 

22 

23 i n s t r 1 

24 kbl = 0.7 

25 kb2 =0.7 

26 ; sinusoidal input 

27 ain oscili 1, 440, —1, 0.25 

28 asig Fbam2 ain,kbl,kb2 

29 asig balance asig , ain 

30 outs asig , asig 

31 endin 

32 

33 </CsInstruments> 

34 

35 <CsScore> 

36 i 1 0 5 

37 </CsScore> 

38 

39 </CsoundSynthesizer> 

For a deeper analysis we can understand the 
system as a 2nd order PLTV system by rewrit¬ 
ing the system equation as 

y(n) = x(n)+/3iai(n)y(n-l)+/3 2 a2(n)y(n-2), 

(19) 

with x{n) = ai(n) = a 2 (n) = cos (toon). This is 
the simplest example but the analysis equations 
are complicated for this second order system. 
However we can also uncouple the input and 


modulation signals as independent streams, and 
treat the whole system as the combination of 
two first order units to reduce the difficulty of 
the analysis. 

Despite the initial results about the 2nd order 
FBAM that were already reported, we plan to 
proceed with more thorough investigations es¬ 
pecially regarding its stability and to determine 
any consequent restrictions it might have on the 
coefficient modulation waveform. 

4 Conclusions 

In this survey we presented some techniques 
for sound synthesis derived from the concept 
of phase distortion. Results can be used to 
generate approximations of classic sawtooth os¬ 
cillators and more timbrally-involved FM-like 
spectra. Implementations in Csound, which 
are easily translated to other languages, were 
presented in order to promote the techniques 
among our community and audio tools develop¬ 
ers. 

The application of time-varying systems is 
not new in general signal processing, but only 
recently is the theory behind these systems be¬ 
ing explored for time-varying digital audio filter 
systems. This paper shows there is a signifi¬ 
cant potential for a number of novel techniques 
based on PLTV systems. The field is still open 
to investigation, in particular with regards to 
2nd and higher order systems. 

We hope that this brief survey invites more 
musicians and technicians to explore the nice 
retro motivated sounds obtained with these sys¬ 
tems, maintaining the interest in distortion syn¬ 
thesis/effects techniques. We hope to encourage 
the sharing of ideas and principles around these 
techniques, such as, for instance, nice distortion 
functions that could be used either in synthesis 
or acoustic instruments processing. 
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